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Abstract 



We evaluate the finite-temperature Euclidean phase-space path integral 
for the generating functional of a scalar field inside a leaky cavity. Provided 
the source is confined to the cavity, one can first of all integrate out the fields 
on the outside to obtain an effective action for the cavity alone. Subsequently, 
one uses an expansion of the cavity field in terms of its quasinormal modes 
(QNMs) — the exact, exponentially damped eigenstates of the classical evolu- 
tion operator, which previously have been shown to be complete for a large 
class of models. Dissipation causes the effective cavity action to be nondi- 
agonal in the QNM basis. The inversion of this action matrix inherent in 
the Gaussian path integral to obtain the generating functional is therefore 
nontrivial, but can be accomplished by invoking a novel QNM sum rule. The 
results are consistent with those obtained previously using canonical quanti- 
zation. 

PACS numbers: 05.30.Ch, 02.30.Mv, 11.10. Wx, 42.60.Da 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Open systems have been the subject of numerous investigations both in classical and 
quantum physics. Physical examples occur for instance in optics — e.g., in cavity QED |IJ 
with obvious applications to laser physics [§], and in microdroplets, which provides the 
spherically symmetric analog — and solid-state physics, where Josephson phenomena || and 
Kondo and related problems @] are all amenable to a "system-bath" description. The open- 
system concept is also relevant to elastic waves and acoustics (e.g., sound emanating from 
musical instruments) and on a very different scale to astrophysics and gravitational waves 

B- 

In a series of papers, we have been particularly interested in those open systems which are 
governed by a wave equation (Eq. fl2.1p below or, by a simple transformation |6j , the Klein- 
Gordon equation [df— d%.+V(x)]ip(x, t) = 0). The "universe" in which the waves propagate is 
supposed to consist of a "cavity" (the open system) and an infinite "outside" , and dissipation 
is caused by waves escaping from the former to the latter J/J. Attention has been focused 
on the role of the cavity resonances or quasinormal modes (QNMs) — eigensolutions of the 
(non-Hermitian) evolution operator of the open system which are exponentially decaying 
in time. It turns out that under certain conditions, to be specified later, the discrete set 
of QNM wavefunctions is complete in the cavity region of space and hence can be used 
for exact eigenfunction expansions. This eliminates the outside from the description; in 
particular, one no longer has to deal with the dense set of modes of the infinite universe. 
In terms of these QNMs, one can establish a formalism which closely parallels the usual 
one for eigenfunction expansions in conservative, Hermitian systems. Applications of this 
formalism include perturbation theory and (canonical) second quantization [§,0, the latter 
being of particular interest to us here. These developments are reviewed in Ref. ||10|| . 

Of course, there are many other ways to eliminate the outside (or bath) and obtain an 
effective description in terms of the damped degrees of freedom of the system only, leading 
to, inter alia, Langevin and master equations PJTl|. A formulation which has proved to be 



especially suited for the study of open quantum systems is the path integral @JT3- First, 
one writes down a functional integral for the generating functional (or density matrix) of the 
universe; the pertinent action, which is readily deduced from the Hamiltonian, is supposedly 
known. Subsequently, one performs the integral over the bath degrees of freedom only, in 
which the system (or cavity) coordinates figure as constant parameters. Since the bath is 
usually taken to be harmonic (in fact, a meaningful separation into system plus bath depends 
on the latter having this or some other simplifying property), this step can be performed 
in closed form. One is left with a path integral over the system degrees of freedom, in 
which the integrand involves an effective action which accounts for the environment in a 
way that is guaranteed to be consistent with quantum mechanics (as opposed to being 
merely phenomenological). This reduced path integral then is a convenient starting point 
for approximations, qualitative considerations or numerical evaluation. 

This paper's aim is, in short, to initiate a synthesis between the QNM and path-integral 
approaches to open systems characterized by escaping waves. To this end, in Section |I| we 
first review the basics of the classical QNM expansion for such wave systems to make the 
paper self-contained, and to establish some notation. In Section [TTT] we present the (finite- 
temperature Euclidean phase-space) path integral for a simple class of wave systems. Since 
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our ultimate interest is in quantities involving the cavity fields, the source function in the 
generating functional is chosen to couple to those fields only, which facilitates the subsequent 
elimination of the bath. This step does not yet involve the QNMs, and the ensuing effective 
action (in particular, the damping term) turns out to have a form which is still well known. 

In Section [IV], the completeness of the QNMs for the models considered is used to write 
this effective action in terms of the QNM basis. This brings together the advantages of an 
effective action with those of a discrete basis [TO]. Up to this stage, the calculation has 
only used the assumption that the outside is harmonic. However, assuming the cavity to be 
harmonic as well, it turns out that substantially more can be done: the special properties of 
the QNMs — in particular, a novel sum rule which we also derive — enable the cavity integral 
to be performed, as an integral over the QNM expansion coefficients. (If the cavity action 
is nonlinear, this integral over the harmonic part is the starting point for the perturbation 
expansion; see Section ]VI|.) Since the effective action is bilinear in the fields and since 
the QNMs are not orthogonal in the usual sense, this involves inversion of a non-diagonal 
infinite matrix (this is to be contrasted with the integral over the outside fields, which can 
be performed for each degree of freedom separately). In our non-interacting problem, this is 
the nontrivial step, brought about by the nontrivial mass density distribution p(x) in (|2.1|) 
below. The result is found to be in agreement with the result of canonical quantization, in 
that both methods yield the same correlation functions. While the damped systems in Refs. 
PIJ^JT^I typically have few degrees of freedom, only their baths being essentially infinite, this 
paper carries out the analogous program of path-integral quantization for a damped field. 

In parallel with a development in the canonical quantization of the system, in Section [V] 
we discuss the special case that the external source couples to the scalar field but not to its 
conjugate momentum. Then, the aforementioned infinite matrices can be transformed into 
a diagonal form, i.e., one involving only a single sum over QNM indices. While the resulting 
expressions look simpler, limitations to their applicability are pointed out. Finally, some 
closing remarks are made in Section W\ 

The appendices contain supplementary material. Since the subjects of Appendices B 
and C are not specific to path integration, they are of independent interest for the study of 
QNMs. 



II. CLASSICAL FIELDS 

In this section, we summarize the QNM expansion for classical fields. In this paper, we 
deal with scalar fields in 1 d only. 

For closed, linear systems, eigenfunction expansions, based on the eigenfunctions or 
normal modes (NMs) of their evolution operators, are a tool of vital importance in theoretical 
physics. However, open systems are not directly amenable to an NM analysis. In these 
systems, any initial state decays in time, so stationary NMs do not exist. As the simplest 
example, we shall be concerned with the real scalar wave equation in one space dimension, 

p{x)dU = d 2 x <t> (2.1) 

studied in a cavity < x < a, with the nodal boundary condition 

0( x =O,t) = O (2.2) 
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at one end but with the outgoing-wave condition (OWC) 



<f)'(a + ,t) = -<j)(a,t) (2.3) 

at the other. The OWC states that, just outside the cavity boundary, the field (fi(x,t) is 
an outgoing wave 0(x — t) (cf. (|2.7| ) below); the condition is specified at a + because, as we 
shall see below, one is often concerned with models in which there is a singularity in p(x) 
at x = a, leading to a possible discontinuity in 0'(x) Jl3|]. The boundary condition ( |2.3| ) 



turns the cavity into a dissipative system that is leaky but not absorptive. The model ( |2.1| ) 
has been widely used as the scalar model of electromagnetism in an optical cavity . More 
physically, the 1-d nature is realized in Fabry-Perot cavities with lengths much smaller than 
the lateral dimensions, and the scalar-field model is rigorously applicable to the transverse 
electric sector. 

For the open system (|2.1|)-(|2~3|), the eigensolutions, labeled by an index j, have the form 

cf>{x,t)=f j {x)e-^ t , (2.4) 
with the QNMs or cavity resonances fj satisfying 

[d 2 x + p(x)uj]}f 1 =0 (2.5) 
and the boundary conditions (|2.2|) , (|2.3| ) translating to 



/ i (0)=0 > f'(a + ) = iu j f j (a). (2.6) 



It is easily verified that Imujj < 0, so that the solution (|2.4j ) is indeed decaying in time. 
Furthermore, the frequencies Uj, which we suppose to be ordered according to increasing 
real parts, are spaced by Au ~ 7r/a, approximately as for a conservative system of size a. 
With the possible exception of modes with Reujj = 0, the QNMs always occur in pairs with 
u_- = —uj*, and one can choose / • = f*. While the field is real, the eigenvalues and 
eigenf unctions are complex; this is the reason for the pairing of modes. 

The usual formalism concerning eigenfunction expansions relies on the hermiticity of the 
evolution operator, which only holds in the conservative case, and therefore breaks down for 
open systems. One possible resolution is to embed the cavity into a universe < x < A 
with a nodal condition at x = A —>■ oo, and study its NMs — the modes of the universe. 
Namely, the system fl2.1| )-( |2~3"D is the restriction to x < a of the problem ( |2.1| ) on the half 
line < x < oo, if one sets 

p(x>a) = 1 (2.7) 

and with the extension of the initial conditions to the "outside" x > a obeying 
(f)'(x>a, t=0) = —<p(x>a,t=0). However, this has the obvious disadvantage of having to 
work with a continuum of states (spaced by Au; ~ n/A — > 0) as opposed to the discrete set 
of eigenfunctions in the conservative case. Besides, the closed system of equations (|2.1|) - ([2.3|) 
shows that even in the presence of dissipation the (thermo-) dynamics of the cavity can be 
studied without explicit reference to the outside, which is the principal goal of the program 
of second quantization of the open system. 
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Previous work (see |10[ and references therein) has established that, in spite of the lack 
of Hermiticity in the conventional sense, an eigenfunction expansion for outgoing waves in 
classical open wave systems can be formulated in terms of the cavity degrees of freedom 
only, overcoming the disadvantages of the modes of the universe approach. The sufficient 
conditions for this QNM expansion are as follows. 

(a) The function p(x) has at least a kink discontinuity at x — a. This demarcates a 
well-defined cavity region from the rest of the universe. In fact, at several places 
the type of discontinuity will be of relevance to us, and we define a parameter p by 
p(x) = p T eg(x) + p6(x — a), where p Te g(%) has at most a step discontinuity at x — a. 

(b) The function p(x) has no tail outside the cavity, i.e., p(x>a) = 1. This condition 
ensures that the outside does not reflect outgoing waves back into the cavity, enabling 
the complete elimination of the environment from the equations of motion. 

These conditions are satisfied for optical cavities bounded from extended vacuum by a 
sharp material interface. Under these conditions, the eigenfunction expansion is exact for 
any amount of dissipation. 

The completeness of the QNMs can be pursued at two levels. First, one shows that the 



retarded Green's function of the system has the representation 14 



G R (x,y;t) 



J2 fj( x )fj(y) e -ta, it 



(2.8) 



for < x, y < a and t > 0, where the /j's are normalized according to ( |2.12|) below. Thus, 
the dynamics is contained entirely in the QNMs. 

Second, realizing that the wave equation ( |2.1|) , like any classical Hamiltonian problem, 
requires both position and momentum to be specified as initial data, one introduces function 
pairs cf) = ((j), 0) T with the conjugate momentum (f> = p<fi, so that for eigenfunctions /, = 



(fj, —ipujjfj) T [T5|]. The set of all function pairs (in general allowed to be complex) satisfying 
the boundary conditions ( |2.2| ) and ( |2.3| ) will be denoted as T — the space of outgoing waves. 



A report on the formal mathematical construction is currently in preparation [jig] 



Using these pairs, one can prove that the time evolution generated by (2^8) can be recast 
in the form 



(2.9) 



where the expansion coefficients are given by 

1 



a j (t) = —(f J ,<f>(t)) 



2uj 



(2.10) 



with a,j(t) = a., (0)e and the bilinear map for £, x G T 

(C x) = i{J dx [CO) xO) + C0*0 xfc)] + C(o) xO) \ ■ 



(2.11) 
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By simply letting t J.0 in ( |2.9|) one arrives at a two-component expansion for an arbitrary real 
(p e T. This expansion makes the completeness of the QNMs manifest. The normalization 
used in ( |2.8| ) to fl2.10| ) can be concisely expressed as 



ifj,fj) = 



(2.12) 



It is seen that ( ]2.12j ) in general is not real, underlining the difference between the form ( [2.11|) 
and a conventional scalar product involving complex conjugation. The fact that (|2.12| ) is 
bilinear also serves to establish a phase convention for the wavef unctions. 
Upon introducing the two-component evolution operator 



Ti 



p{x) 
dl 



-1 







(2.13) 



the cavity evolution ( p.l|) can be written as id t <p = 7i,<p, in striking analogy with quantum 
mechanics. In this notation, the definition ( |2.5| ) of fj takes the form Tifj = ojjfj. The 
operator 7i can be shown to be symmetric with respect to the form ( |2.11| ), i.e., 



(2.14) 



for any £, x £ T. This analog of hermiticity holds even though the system is not conservative. 
The symmetry of Ti yields the "orthogonality" relation 



(fj, fk) = for ujj ^ u k 



(2.15) 



in an immediate transcription of the usual proof, leading to the uniqueness of the expansion. 
Incidentally, an expansion such as ( |2.9| ) but involving the first component alone would not 
be unique. 

Instead of its present formulation as an "orthogonal" expansion involving a bilinear map, 
the series (|2.9| ) can also be regarded as a bi-orthogonal expansion in terms of the standard 
inner product |T7|]. The power of this latter, slightly more involved method shows when 



several QNMs merge | lq| , a possibility which in this paper we will only briefly consider in 
Appendix A. 



III. ELIMINATION OF THE OUTSIDE 



We want to calculate the generating functional for the cavity field, in a form which 
manifestly involves the QNMs. Since we want results for finite temperature, a Euclidean 
formulation is advantageous. The QNM expansion ( |2.9j ) involves two components; therefore, 
we have to use the phase-space version of the path integral: 



S{X} = (r T ex V yyr(cf>(T),x(r))}\ 

= Z- 1 Jv ( (>(x,T)expUjT (<I>,X)-J Q 



/(, \2p 2 



(3.1) 
(3.2) 



In (|3.1|) , x is a real external source ||19|| ; the coupling to the cavity field only (cf. the In- 
troduction) is taken to be of the form ( [2.11|) , anticipating that, by (|2.15|) , this will simplify 
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expressions upon QNM expansion as in ( |2.9| ). The inverse temperature is (3 = 1/T (through- 
out % = k-Q = 1). Imaginary-time ordering is denoted by T T as usual; this operation is needed 
in (|3.1| ), in which <p is still a Heisenberg operator field, but has no place in ( |3.2| ) which in- 
volves c-numbers only. In the following, the meaning of <f> will be clear from the context. 
The two-component measure is V<p = X?0T>0. The normalization Z equals the path integral 
on the r.h.s. with x l— *• 0- This is to be understood formally, for in fact both the integral and 
Z are infinite, only their ratio being meaningful. In the course of the calculation we shall 
cancel all factors which do not depend on x against corresponding factors in Z, without 
bothering to reflect these subsequent redefinitions of Z in the notation. 

The boundary conditions on the real field <j> are 0(0, r) = 0(a + A, r) = 0, <p(x, 0) = 
4>(x, /3), the latter arising from the trace implicit in the quantum statistical expectation value 
( p.l| ). No boundary conditions can be imposed on 0, since the absence of terms with 0' in the 
action implies that typically is completely discontinuous. This at the same time means 

J2Tf. However, these problems in 



that phase-space path integrals can be quite tricky 
general show up only beyond the semi-classical approximation, and hence should be absent 
here since we look at a linear problem for which this approximation is exact |[22| . We are 
interested in the limit A — ► oo, but postpone this operation until it can be performed in a 
controlled way. 

We want to split the integral into a cavity and a bath factor: Z~ l ^Dcfi = Z~ 1 j'U(f) c x 
Z^ 1 JT><p h , where the latter runs over fields on (a, o+A), with a given boundary value 0(a, r). 
The integral over bath momenta is trivial; introducing £ = x — a, 



V<t>h 



exp 



dr f d£ ( — — h — i 

o io S 2 2 




exp 




(3.3) 



where we used fl2.7p . The 0t,-integral is nontrivial only because of the boundary condition, 
which can be implemented through the expansion 




sin 



TTU^ 



(3.4) 



with the bosonic Matsubara frequencies v m = 2nmT, m G Z. 
Substituting (|3.4j) into (|3.3|), the action can be written as 




where 



A" 1 |0 m (a)| 2 + E 

/ u=l 

\ oo 

A- 1 |0 m (a)| 2 + £ 



u=l 



+ 



2A 



2 2A 



+ R e (0 um 0_ m (a)) 

Till 

2(A^/™) 2 |0 m (a)| 2 



Az/^ + 7T 2 M 2 /A 



(3.5) 



+ 



2vl l cj) m {a)/'Ku 

V 2 m + (7I"M/A) 2 



(3.6) 
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Switching to variables cj) urn does not affect the domain of integration (in particular not in a 
m (a)-dependent way), since both the (p um and the <fi um run over all C, subject only to the 
restriction (p um = 0* _ m and analogously for 4> um . The completion of the square in the last 
line of (|3.5| ) thus eliminates the need to actually perform the path integral, and one obtains 
some constant independent of 0(a), which then cancels against the same factor in Z^. This 
leads to 



V(f) h 




A 



-i 



2(A^/7m) 2 

^1 A^ + 7T% 2 /A 



E 



exp|-^^|z/ m ||0 m (a)| 2 | for A 



oo, 



(3.7) 
(3.8) 



where to arrive at the last line we evaluated the sum over u asymptotically in A 1 , using 

1 



5 u 2 (l + e 2 u 2 " 



oo 

E 

u=l 


L 2 


7T 2 




6 


-/ 


7T 2 


7T 


6 


c 



u 2 + e 2 
du 



o m 2 + e 



+ 0(e 2 



(3.9) 



for e = 7r/Az/ m , leading to the cancellation of the O(A) terms in the exponent of (|3.7|). 

The Caldeira-Leggett type 0] exponent in (|3.8|) is the quantum-mechanical, finite- 
temperature equivalent of an Ohmic-damping term. Its emergence here is not surprising, 
given the correspondence between transmission-line environments of the type considered 
here and the oscillator baths used in its original derivation [FJ|. Namely, in the limit 
A — > oo waves escaping into the homogeneous outside string will never be scattered back, so 
that the outside acts as a sink for waves emanating from the cavity. Since the model ( |2 . 1|) 
is dispersionless, this damping is frequency independent. On the classical level this leads to 
( p.3|) , in which on the l.h.s. is precisely the string tension; according to the r.h.s., this force 
is equal to —0. This velocity-proportionality is reflected by the first power of v m in (^.8|); 
however, unlike (|2.3|) , the action ( |3.8|) is necessarily (imaginary-)time reversal invariant. 

Substituting the bath contribution back into the generating functional ( p.2| ) and trans- 
forming to Bose frequencies also in the cavity, one gets 



SM = /^ expiry 



Wn 



(,4>m,X-m) 2~l0m(a)l' 




mfrn'r —m 



2p 2 

This form completes the elimination procedure in that it manifestly involves <p c only. 



(3.10) 
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IV. PERFORMING THE CAVITY-FIELD INTEGRAL IN THE QNM BASIS 



A. Path integral and Green's function 



To make further progress, we introduce the QNM expansions 4> m = J2j a jmfj an d Xr, 



Y^,jbj m fj f24|1 ; the coefficients satisfy 
and invoking the "orthogonality" relation 



b*_j _ m . Substitution into ([3.1 0|) 



/ dxpfjf k = 5 jk -i 

•/ 



fj{ a )fk{o) 
UJj + UJ k 



(4.1) 



which follows from ( |2.11| ) and ( |2.15| ), leads to 



S{x} = /^exp{TE 



E 0'jmO'k-m,fj{.o)fk{o) 
jk 



v m (9(m)uj k - 8(-m)ujj) + iujjUJ k 
UJj + uj k 



+ E 2 ^ b 



j,—m^jm 



(4.2) 



ex PS T E 



2J ajm.ak-mfji^fkiO' 
jk 



v m (9(m)uj k - 9{—m)u)j) + iujjuj k 



~ E b imk 



UJj + UJ k 

fj(a)fk(a) ( 9(m)uj k 9(-m)ujj 

k,—m , I • . i . 

UJj + UJ k \lUJ k + V m lUJj - v m 



(4.3) 



where 



Ijm + E 



km ' 



fj(a)fk(a) ( 9{m) 



ujj+uj k \uj k (iujj + v mj 



+ 



9{-m) 



UJj (lUJ k - v m ) 



(4.4) 



and 8(0) = ~. To arrive at ( f4.2|) , one has to use the fact that the Kronecker term in ( |4.1|) 
does not contribute in the action by m-parity. Thus, it is precisely the second term in ( f4.1| ) 
which contributes; its presence makes fl4.2| ) non-diagonal (i.e., involving a double sum J2jk)l 
the surface values fj(a)f k (a) are a measure of dissipation, since they would vanish if the 
system satisfied a nodal boundary condition also at x = a J25 ] . 

The only nontrivial ingredient in the completion of the square (|4.3| ) is the QNM sum rule 



E 



fj(a)fi(a)f e (a) 

k (UJj + uJ k ){uJ k + UJl) 

which will be discussed in Section II V B. 
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(4.5) 



26 



We thus arrive at the final answer for the generating functional in terms of the coefficients 

fj(a)f k (a) 



S{x} = expl %T V b jm b k ,„- 



9(m) 



u k 



+ 9(-m) 



UJi 



UJj + iv m t 



(4.6) 



^k — %v m 

The relation between S{x} and the temperature Green's function Qj k {r) = — (T T {aj(r)a k }) 



is standard [27 
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d 2 S{ X } 



dbim.dk 



'jmUUkj—m 



db jm db kr 



d 2 I r rP 

T T expiry 2wjb jm j dr e~ lVmT t 

\ jm 

Aw^T 2 l P dr x dr 2 e^-n) (^(nKfa)}) 

J 

-AujU k TQ jk {v_ m ) , (4.7) 

that differentiation (in our case 



where it is a standard result of functional integration 
ordinary partial differentiations with respect to the discrete basis {bj m }) of the generating 
functional automatically yields time-ordered expectation values (cf. below ( |3.2| )). Substitut- 
ing (|4.6|) into ( |4.7|) , one obtains 



—i 



fj(a)fk(a) 
2u} j u} k {u j + to k ) 



9{m) 



LUj — lV n 



6{-m) 



u k 



u) k + iv n 



(4.8) 



Quite generally the temperature Green's function is related to the real-time retarded prop- 
agator Gf k (t) = -i9{t)([aj[t),a h ]) by 



G%{iv m ) 



(4.9) 



for m > 1 f29f . Evaluating the analytically continued G R k at the discrete frequencies 



( |4 . 8| ) thus is readily seen to be in exact agreement with the results obtained in Ref. |J by 
canonical quantization. 



B. Action-inverting sum rule 

In the previous subsection we deferred the derivation of (|4.5|) not to disrupt the flow 
of the argument; it will be discussed presently. For j ^ £ one can apply a partial- fraction 
expansion to the summand and see that the sum indeed vanishes if J2k f k { a ) / i u k + Wj) is 
independent of j. This can be shown by rewriting 



j 



h UJ k + LUj ~ \LU k U k (iO k + U)j) J 

j2 $W + 2^G R (a, a; 

k U k 

+ (4.10) 

k U k 



where the last line follows because more generally one has 



2uo J G R {x ) a ] -uo 3 )=i f -j^\. (4.11) 

For a proof, let to — > — uij in the purely classical identity G R (x, y\ uj) — G R (x,y; —uj) = 
(2uj/i)G R (x,a;uj)G R (y,a] —uj) established in Ref. ||, and compare residues on both sides. 
The value of the first term on the r.h.s. of (|4.10|) is irrelevant for the derivation of ( |4.5|) ; in 
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Appendix B it will be shown that, for the case of a step discontinuity in p(x=a), one has 
E k f 2 k (a)/(u k + uj) = i(p(a-) + l)/(p(o") - 1). 

It remains to prove ( |4.5| ) for j = £, for which one has to calculate — J2k fk( a ) I i^k+^j) 2 = 
2d u) [ujG K {a,a]oj)} UJ= ^ U j.. To this end, define f(x,u) (g(y,u>)) as the solution of the homoge- 
neous wave equation ( |2.5| ) (upon the substitution Uj \— > uj) satisfying the first (second) of 
the boundary conditions (|2.6| ). This allows one to write 

CV^)= ^^ , (4.12) 

where one can choose 

f(x,u)=f(x,-u), (4.13) 

and where W = fg' — gf is the position-independent Wronskian of the functions / and g 
\T§. Together with the OWC for f(x,Uj) = fj(x) and g, ([112]) yields 



2uj j f j {a) 



If Uj were a double QNM pole []Iq1 , i.e., if / satisfied the OWC up to and including 0(uj—LUj), 
the numerator would vanish. However, we assume simple poles throughout the main text 
(cf. Appendix A), so further evaluation is needed. Solve [d% + pu 2 ]d u} f\ u}j = —2ujjpf. 
by the variation-of-constant Ansatz d^fi^x.ujj) = fj(x)hj(x), leading to h'-(x)f 2 (x 
2u j fZdyp{y)ff(y) h'{a+) ^ i-(fj, ^/fKa) ^ iu j d^f(a,u j )+if j {a)-d u f(a+,u 



3 j 



(fj, fj)/fj{a); given our observation that (^4.14|) vanishes for a double pole one could have ex- 
pected an answer oc (fj, fj) |Tj|, which in the present case of course equals 2uj. Substituting 
back into ( |4.14j ), one finds 

2d^uG K (a,a;u)] u= ^ 3 = (4.15) 

Jj \ a ) 

completing the proof of the sum rule ( [4.5|) and therefore of ([Of). 

V. ONE-COMPONENT FORMS 

Given previous experience with the canonical approach 0, one would expect that the 
cavity functional and effective action also have a diagonal form (cf. Appendix C) if the source 
couples only to the first field component <fi. In this section we thus temporarily set the first 
source component to zero: x — 0, and consider S{x}- For x = 0, the sum rule (f4~5p implies 
that the bj m satisfy 

\ - fj(a)fk(a) , (KU 

2^ 1 b km = -lOjm , (5-1) 

k Uj + Uk 

which at once yields 
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S{x} = ex p| T E & i 



h ■ 



jm 



UOj — 1 1 is n 



(5.2) 



However, this form should be used with caution: the fj are overcomplete for the one- 
component expansion of x, and consequently the bj m are now no longer independent, as 
exemplified by (|5.1|). It might thus be more convenient to write out the projection formula 



V 2a; j)Jo a dxfjXm, i-e. 



S{x} = ew\ -7E / dxd VXr 

{ 4 im J 



{x)x- m {y) 



LOj(ujj - i\v m \) 



[ T 



[ l m Jo 



(5.3) 
(5.4) 



Eq. ( |5.4| ) looks very simple, yet this form (in particular the representation (|5.3| ) of Q) would 



have been difficult to derive from (|3.10|) without using the power of the two-component 
QNM expansion in the intermediate steps, cf. below ( |5.8|) . Using the representation ( |5.3j ) for 
S{x} = (T T exp^i Jq(1t Jq + dx x(x, r)<p(x, r )}), 4>~4> correlators can be found by functional 
differentiation with respect to x, without the problems associated with ([5.21). 

One can also obtain a diagonal representation for the effective cavity action. This might 
for instance be useful in the numerical study of an extended, interacting, model; cf. the 
discussion in Ref. ||, Sections VI and VII. For x — 0, one can integrate out and study 
the configuration-space path integral 



/ z exp l 


m 


[We J 

/ z exp 1 


m 



| 2 -^IC! 2 



j jk 



2{uj + u k ) 



-fj( a )fk(a) 



The aj m are now of course given by the projection 



(5.5) 





\ r 





2cu? 



+ 



(5.6) 



implying 



x - uj k fj(a)f k (a) 

2 , : a kr. 

k 10 j + LOk 



tUJjCLj m 



(5.7) 



By also using J2j fj( a )/ u j = (see ([Bo])) and the representation m (a)0_ m (a) 
— i(4> m , <fi~ m ) = —2iYlj ^jCtjrn a >j,-m, the one-component action can be written as 



12 



s{x} 



v<t> c 



exp \ T ' 



l jra 



l\V m \UJj — iav m ujj 



(5.8) 



j'" 



Like for the one- component functional Q5.2D , one has the caveat that the a jm now no longer 
are independent degrees of freedom (as, apart from the constraint a- m = a*_j_ m , they were 

is difficult, and the best route to (15. 21) 



in Section [TV|) . Therefore, direct evaluation of 
is probably via the auxiliary field as has been done in the preceding. Again it might be 
more convenient to write out the expansion coefficients, viz., 



s{x} 



dx i 



iX—m 



4 Jo 



dxdy <P' m (x)<p'_ m (y) 



x 



E 



vt + ui — i\v m \ujA - iavLu 



mm 



(5.9) 



VI. DISCUSSION 

By carrying out the path- integral quantization of the open system Q2.1| )-( pT3"D , we have 
met the challenge set out at the end of Ref. [[|. This leaves us facing the second chal- 
lenge mentioned there: the inclusion of matter in the system's Hamiltonian or action |3~0 



Assuming the interaction to be confined to the cavity, the elimination procedure of Sec- 
tion H| goes through unmodified and one simply has an extra term in the effective-action 
exponent of ( |3.10| ). It has already been mentioned in the Introduction that the analysis 



of Section [IV A| consists of two steps. The first of these, expansion of the effective cavity 



action with respect to the QNM basis, certainly goes through for any form in the interac- 
tion since the QNMs are complete. One obtains an appropriate generalization of ( f4.2|) , and 
the advantages of a discrete basis now apply to any qualitative, approximate or numerical 
analysis. While the second step of closed-form evaluation will not be possible any more in 
general, by substituting dj m i— > d / d{2ujjTbj^ m ) in the interaction term one has the stan- 
dard representation of the interacting S{x} as a functional of the free 5*o{x} [p8| , the latter 
being given by (|4.6|) . For instance, if one incorporates into the cavity action a toy nonlinear 
term S mt = Jq dr J a dx A(x)0 4 (x, r) (so that the integrand of (|3~^ ) acquires an extra factor 
exp{— <Si nt }), one can write 

S{ X } = Z- 1 exp{-/3 £ 8 mi+ ... +m4 £ Xj 1 ...j 4 If oTST-iw , (6-1) 

I mi...m 4 ji-.-ji i=l 3i O0 ji™i ) 

with the coupling A^...^ = J a+ dx\(x)fj 1 (x) ■ ■ ■ fj 4 (x). This formal representation is a suit- 
able starting point for perturbation theory. 

It should be noted with some modesty that the final result ( f4.6|) could have been known 
in advance even without employing results of the canonical quantization of the system (cf. 
0, Section III). Namely, the kernel (|4.8|) is the analytic continuation of Gf k (u) by ( f4.9|) . In 
their turn, the Gj- k (t) are the unique QNM expansion coefficients of the tensor correlator 
G R (x, y; t) = —i9(t) ([4>{x, t), <f)(y)}), all components of which follow by trivial time differen- 



tiations from its (1, 1) component (cf. ( |A24j) ), the latter merely being the classical Green's 
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function ( |2.8|) ; see (||, Section VI B and Appendix C) pi] . Yet, the present explicit eval- 
uation of the path integral is of considerable interest, since so few of them can be done in 
closed form unless they trivially decompose into a product of ordinary integrals over normal 
modes. Also, carrying out the calculation has uncovered new results on QNMs — the identity 
([4.5|) and those of Appendix B — which are useful already on the classical level. 
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APPENDIX A: JORDAN-BLOCK PATH INTEGRAL 



1. Evaluation of the integral 

The results in the main text have been derived under the assumption that all poles in 
G R (u) are simple, as is for instance obvious by considering the Fourier transform of 



Here we consider the general case, and follow Ref. [18| throughout ||32|| . For each QNM pole 
in G r (lu) of order Mj at uj = ujj, one introduces 



fj,n{ x ) 



(Al) 



for < n < Mj — 1, with f(x,u) defined above ( 4. 12| ) . The associated momenta read 

fj,n = ~ip[uJjfj,n + fj,n-l] (A2) 

(fj n = for n < —1), so that fj t0 = fj is the QNM eigenvector, which together with 
{fj,n}n=i 1 spans a so-called Jordan block of the Hamiltonian (|2.13|) . 



One of the main results of 181 can now be stated as follows. Provided that one chooses 



the functions / and g such that, for all j, 

g(x,u) = /(x,^) + (9[(^-^) Mj ] 

and 



(A3) 



(A4) 



(where we draw attention to the orders of the error terms), which is readily achieved, one 
has the biorthogonality relation 



\fj,ni fk,r) 

for < n < Mj - 1, < r < M k - 1. 



2ujjSj k 5, 



n+r,Mj-l 



(A5) 
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j Z^n=0 u j,n;mJ j,n 



S ubstit uting the expansions <p m = Ej En=o aj,n;mfj,n and Xm = Ej E 
into ( |3.10|) , the term involving /' n can again be dealt with using an integration by parts, 

supplemented by the QWC f^ n (a + ) + fj,n(a + ) = and the wave equation [d% + pujf\fj^ n = 

—p[2ujjfj, n -i + fj, n -2\- However, integrals J a+ dx pfj, n fk,r m general cannot be reduced to 
surface terms by a single application of ( |A5|) as in ( |4.1|) , because of the second term in ( |A2|) . 
Rather, iteration leads to 



5 jk 6(n + r + I - MA " ' 

dx pfj, n f k , r = , ,, Wr4 ,1_ M ., + * E E 



-U j ) n + r+1 - M i 



p=0q=0 



P + q \ fj,n-p(a)fk, r-q(a) 
P I {-LOj - UJ k 



\p+q+l 



5 jk 6(n + r + 



Mi 



32 k 5 £k f{a,u)f{a,p) 



-u 



\n+r+l-Mj 



n! r! 



Co> + /! 



(A6) 
(A7) 



where the equivalence of the second line to the first follows from (|A1| ). Comparatively 
compact forms like ( |A7| ) will be essential in the calculation below. Comparing the second 
(non-diagonal) terms on the respective r.h.s.'s of (|4.1|) and fl^7f) , the latter is seen to be a 
differentiated version of the former. However, it should be borne in mind that ( |A7|) only 
holds for the special choice of / satisfying ( |A3|) and (|A4|) . 
Returning to ( |3.10| ), the above yields 



S{X} 



Vcp c 



^ ] ^ ] ^ ] ^'j,n;m^'k,r;—m^j,n;k,r(j^m) 
jk n=0 r=0 



+ E E 2w 

j a j,n;mbj,Mj-l- 

j n=0 



with the action matrix 



< Sj,n;k,r^m) 



an I I 

|Wj ^ \U>k 



n\ r\ 



f(a,u)f(a,p) 



v m (pQ{jn) — u9(—m)) + iup 
io + p 



It is seen that the first term of 



cancels in S, cf. the remark below 



We claim that the result of (IA8I) reads 



(A8) 



(A9) 



S{x} = exp<{ -2T E E ^j^kbj^mbk^-mGj, 

jkm n=0 r=0 



in terms of the Jordan-block temperature Green's function 

Q j,n\k,r{ym 



dMj — l—n i 
A I 



(M,-l-n)!(M fc 



^ f(a,X)f(a,fi) \ 6{m)\ | 6(-m)p 



r)\2u)jU) k (\ + p) liA + z/„ 



2yU — z/„ 



To verify this claim, it suffices to show that |26| 

1 



M k -1 



M fc -1 



E E <WW£W )r K 



SjiSnu — E E ^A,r;j",n( l/ m)^Jfc,r;^,«( 



fc r=0 



fc r=0 



(A10) 



(All) 



(A12) 
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where in view of the simple symmetry Sj, n -k, T (y m ) = Sk, r ;j,n(—v m ) and likewise for Q, the 
second of these relations is equivalent to the first. Therefore, let us evaluate the l.h.s. of 
( |A12| ), first assuming v m > 0. Substituting ( |A9| ) and ( |A11[ ), the sum over r in ( |A12| ) can be 
performed using the product rule for differentiation, and one obtains 



A4-1 _ 
k r=0 







QM t - 



l-u\ 



v m + ioj X 



n\ (Mi — 1 — u)\ v m + iX uj£ 
Carrying out the partial-fraction expansion 

H 1 



f(a,u)f(a,X)^2 



Aft. -II 



\t»>h 



(fi + u)(fi + X) uj - X 



U 



(M k -l)\%) k ((i + w)(ji + \) 



X 



(A13) 



fl + U [L + A 



(A14) 



one recognizes in (|A13|) the QNM expansion of the Green's function Jl8 

I-'* f(x,/j>)f(y,ti) 



QM k -X\ 

g r (x, y; = E {Mk -iy. 2u k (c-^ 



(A15) 



One can now replace (|A15|) by the alternate form ( |4.12|) (which holds independently of the 
pole structure of G R ), and in the latter use ( |4.13| ) to write 



G R (x, a; —u) 



The denominator can be simplified by 



f(x,u) 



f'(a+,u) + iuf(a,u) 



(A16) 



iuf{a, u) - f'(a + , u) = 2uJj{uJ ^ ^ 3 + 0[(u - Uj) 2M * 



g{a,u) 



2u)j(u — Uj) 



f(a,u) 

where we subsequently used ( |A4|) and (JA3|), leading to 

if(x,u) Ujf(x,u)(u-Uj) M i 



+ 0[(u-Uj) 



2M, 



uG (X, a: —uj) 



2f(a,u) 



+ 



2uf 3 (a,u) 



+ 0[(u- Uj 



,2M,- 



(A17) 



(A18) 



which is seen to generalize both ( JOT ) and ( |4. 15| ) simultaneously. The last (higher-order) 
term of ( [A18Q is readily verified not to contribute in the final result ( |A20|) below (however, 
for this it is essential that it be at least of the order indicated), and from now on will not 
be written explicitly. Substitute ( |A18| ) (with x i— ► a, and with (u,j) i— > (X,£) in the second 
term resulting from (|A14|) ) into the upshot of (|A13|) , one arrives at 



E E WWWW = 

k r=0 

df i ~ 1 ~~ u \ ult u m + iuj X f(a, u)f(a, X) 



d n \ 



n\ (Mi — 1 — it) ! v m + iX ui uj — X 



u e (X-u e ) M * Uj(u-Uj) M * 



2XP(a,X) 



2uf 2 (a, u) 



. (A19) 
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1 and u > 0, one immediately sees that this vanishes if j ^ £, in which case 



Since n < Mj 

the denominator uo — A does not become singular. If j = £, in the first term of ( A19 ) one can 
use {y m + iu;)f(a,u))/[(u m + i\)f(a, A)] = 1 + 0{uj — A), where the higher-order term does 



not contribute because it cancels the singular denominator u — A, upon which (A 



yields zero in the final differentiation. The second term of ( |A19| ) is handled analogously, and 
one is left with 



Af fc -1 



E E <WK)£W>„ 
k r=0 



(X — U)j) M i — (u — U)j) M i 



n\ (Mj-l-u)\ 



2(u - A) 



l^k^A k_ £ (cj - UJj) p (X - 0J j ) Mj ~ 1 ~ P 



2 n \ (Mj-l-u)\ £ 



— ~Sn 



(A20) 



which proves the claim made above ( |A10| ) for the case v m > 0. The < is fully 

analogous but slightly simpler, since the factors involving v m cancel from the outset. 



2. Comparison with canonical quantization 

Let us finally make a brief comparison with the operator version of the quantum theory, 
for an independent check on the above algebra and for a better understanding of the unusual 
quantum excitations corresponding to the higher-order poles in the correlation functions. It 
turns out that the generalization of Refs. |8||| is straightforward. Expanding the operator 
field (f){t) = Dn=o 1 a j,n{t)fj, n , the operator QNM expansion coefficients a 3 -, n satisfy the 
coupled (for Mj > 2) system of quantum Langevin equations 

i 

a j>n + iujdj^ + ia jtn+1 = - — f jtM .^ n (a)b(t) , (A21) 

for < n < Mj — 1 (with aj t M 3 = 0), with the environmental driving force b(t) = 20 in (a+t) = 
(f)'(a + ,t) +0(a + ,t) (Eq. (4.2) in M), which is independent of the QNM pole structure. 
The hierarchy ( |A21|) yields dj^^) — (i/^j)fj,M -i-n( a ,^)b( u} ) i n terms of the func- 
tions f\ n (x,u) = f^dte^% n (x,t) = j °°dt e^(n\)^d;\ u .[f(x, fi)e~^] = i ££ =0 fj,n- P (x) x 
(a; — ujjY^ 1 = (i/ntyd^lujlf (x, fi) / (lu — /i)] (cf. |TJ] for fj ;n (x,t)). For the correlators this 
implies 

Fj t n;k,r{ U ) = \%,n( a ') a k,r) 

_ /,, tf; -i-„(«,4,M r i -j^ {l{u)bh (A22) 



where the driving-force correlator reads (b(u)b) = 2u/(l — e' 13 ^) as for simple poles 
(0, Eq. (6.3)). 

The retarded propagator may be found from ( |A22j ) as G^ n;A , r (c<;) = $(duj' /2n) x 
(u-u/ + ley^Fj^iuj') - F k ^ n (-oj')} (cf. §, Eq. (7.2)), leading to 



17 



, \ i d x J 1 n l 3 1 % Xf(a,X)f(a,fi) 

^j,n;kA u ) 2uju k (Mj - 1 - n) ! (M k - 1 - r) ! (u - A) (A + //) ' 1 J 

These are the QNM expansion coefficients of the tensor correlator G R (a;, y;t) = —id(t) x 

<£(?/)]) = EjfcE„io" 1 ErlV lG '^„;fe,r( t )/?>( a; )/A:,r(l/)- Since the latter is related to 
the classical Green's function by 

e 'W>= u$&) ew °- (S **>*(*-»>) (A24) 

(cf. Q, Eq. (6.6); the second term in ( |A24j ) comes from time differentiation of the step 
function in the definition of G R ), (|A23 ) can also be obtained directly from the classical theory 



by substituting ( |A15|) into ( |A24j ) and subsequently applying the tensor-product projection 



formula (i.e., Appendix C of 0, generalized to the present case of nontrivial Jordan blocks, 
which merely involves using ( |A2|) for the second components and keeping track of the intra- 
block indices as stipulated by (|A5|) ). Conversely, using the sum rules implicit in ( |A18| ), the 



expression G R (x, y; u) = J2jk E*=o 1 Er=o 1 G ^n ; fc,rM/i,n(^)./fc,r(?/) is readily reduced to the 



form ([AT5D . 



Since inspection shows that Q as in (|A11|) and G R as in ( |A23| ) are related as in ([D| 



the path integral and canonical approaches to the quantum theory agree for arbitrary QNM 
pole configurations. 



APPENDIX B: QNM SUM RULES 

In this appendix we give an overview of some sum rules which have been used in the 
above, as well as more direct proofs of some others which can be read off by comparing the 
various forms for S{x} in the main text. A central role is played by the WKB approximation 
to the Green's function Ijlif, 



sm(uT(x)) e~ iujT ^ + R(u)e iujT ^ 

G R {x <y<a;u)^ 1 — } - , (Bl) 

uj^Jn(x)n(y) [e- iujT + R(oo)e iujT ] 

where n = y/~p, and with the geometric-optics transit time T(x,y) = J^dsn(s), T(x) = 
T(0, x), T = T(a). For x < a < y one has (exactly) G R {x, y; u) = G R {x, a; uj)e iuJ ^~ a l If the 
discontinuity at x = a is at least a step, the reflection amplitude is 

R(u) = — ; B2 

for weaker discontinuities we write R(u) ~ R p /u p for some integer p > 1. 

Standard wisdom states that G R (t=0) = 0, but in fact this is only true distributionally, 
and pointwise for x^y. Examining the spreading-plateau solution for G R (x,y;t) for small 
positive t (where p(x~y) can be taken constant), one finds that G R (x, x;0 + ) = —l/2n(x) 
if n is continuous at the point x\ at x = a one finds instead G R (a,a;0 + ) = —\/{n~ + 1) 
(n~ = n(a~)) or G R (a,a;0 + ) = if fi > (other possible points of discontinuity inside 
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the cavity work analogously). Since certainly G R (t=0~) = 0, by the general theory of 
Fourier integrals one expects G R (t=0) to be the average of these values, i.e., G R (x, x;0) = 
—l/2(n(x~) + n(x + )); this can be verified by performing the countour integration of G R 
in the upper half plane, where only the large semicircle contributes. Subsequently, one can 
integrate in the lower half plane, which upon comparison yields the sum over QNM residues. 
Define 



s{x,y-t) = Y, f ^ X)f ^ V) e-^ t 



10 A 



(B3) 



since the terms of this sum only tend to zero if t > T(x) + T(y) — 2T/p M, and since we 
only want to consider (|B3|) for t = and t = + but for all x, y in the cavity, the analysis 
will be restricted to p < 1. If t — 0, (|B3|) in general has to be understood as a principal 



value lim 



E 



M 



-M- 



Choosing x < y, one has the following cases. 



x = : 

< x < y < a : 
0<x<a<y<a + T(x, a) 

< x = y < a : 



(all terms zero) 



n[x) 



x 



y = a, t = : 



0. 



2i 



(n 



= 
= 0+ 

li > 
, step 



x 



a, t 



) 2 -l 
divergent, p = 1 

0, ft > 

step 

1 



2% 

'n- + V 



(B4) 



V 



0<x<a, y = a + T(x, a), t = : 
0<x<a, y = a + T(x, a), t = + : 
0<s<a, y = a + T(x, a) + e (e j 0) 



0. 



n 



n — 1 V n{x) 
t divergent, 







0. 



2% 



n 



n — 1 V n^xj 
divergent , 



fi > 
step 
p = 1 

/i > 
step 

p = 1 



Some comments on this table: note that results for /i > follow from those for a step by 
formally setting n~ — > oo, and those for p = 1 follow by letting — >• 1. The result for 
x = y = a, t = + , p = 1 also follows from the one for < x = y < a, t = + , p = 1 by 



letting x ] a (since n = 1 if p = 1). In this latter case, the factor e 



serves not only 



as an oscillating regulator, but also as a power-law damping. If there is at least a step at 
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x = a, however, the sums with t = converge faster than the regulated ones, since in the 
absence of a factor e _lWji the leading 0(j~ l ) tails in the summand for positive and negative 
j cancel exactly. These results can be checked for the "dielectric rod" p(x<a) = const ||, 
where they become conventional Fourier series. 

With (IFJ3D at our disposal, we can return to ( (4.10| ). In the second term on the r.h.s., the 



summand has a denominator ~ lu^, so that one can indeed take a; | a in (|4.11|) , as has been 
done in the main text. In the case of a step in p(x) at x = a 
( |4.10|) quoted below ( |4.11| ); for p > 0, one has J2k fk( a ) I '(^k + u j 



leads to the result for 
= i instead. For p = 1 

one finds that J2k fki^/i^k + u j) diverges, but this merely means that the partial- fraction 
expansion applied to ( f4.5| ) in that case is valid only term by term. In fact, one can readily 
verify that (|4.5| ) itself remains convergent and valid also if p = 1. 

-2G K (x, y; 0) for < x < a, x < y < 
y. For uo = 0, the differential equation for G R is trivial, and one finds 

fj( x )fj(y) 



A different QNM sum is £j f j (x)f j (y) ju) 



a + T(x, a) and x 



E 



2 mm(x, y) 



(B5) 



Operating on this with d\ reproduces the familiar p(x)J2jfj( x )fj(y) 
< x, y < a (more distributional sums are given below). 



25(x — y) for 



The form ( |5.9| ) for S{x} also prompts one to study J2j fji^fjiv)/^ an d J2j /■ 



x) x 



fj(y)/Uj] I restrict myself to < x, y < a. In both cases one can operate with d\ and 
dy termwise for < x, y < a, and determine these derivatives from (|B4j ) and (p5|). Since 
furthermore the sums vanish for x = or y = 0, one has 

C 3 xy 



E 

j 


f:, 








E 

j 


fi 


[x)fj(y) 







C^xy + 2 / ds p(s) min(s, x) min(s, y) 
Jo 



(B6) 



The constants follow from 



d x E fj ^ fj ^ 



3 



C 3 y 



• J2 li ( a )fj(y) , ^ fj ( Q ) fj (y) 



3 



= 2iy 
C, = 2i , 



(B7) 



and 



0»£ 



fj(x)fj(y) 



c, 



(C 4 + 2p)y 

fj( a )fj{y) 



E 



2(n-a)y 
-2a . 



+ /^E 



fj(a)fj(y) 



(B8) 
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These classical sum rules combine nicely to yield some distributional identities. We 
already know p(x) J2j fji x )fj(y) — ^( x ~y) ii < x,y < a, but what if y | a? This question 
often arises in the course of calculations. However, we know little about QNM sums outside 
[0, a], and we want distributions to operate on functions defined in that interval only. Thus, 
we want to evaluate q(x) = J2j 9(a + — x )p( x )fj( x )fj( a )', n °te how writing the discontinuous 
(singular for p > 0) 9p after the sum avoids the need to multiply this factor with the 
singular (discontinuous for p > 0) J2j fj{ x )fj{ a )i so that as written q(x) is well-defined. 
One evaluates 



q(x) 



j 

3 

+ d x ^2(d x 9(a - x)) 



2fj( x )fj( a ) 



U)1 



fj(x)fj(d) 



UJ 



+ ^(46»(a + - x))d x 



f j (x)f j (a) 



J j 
f j (x)f j (a) 



^'7 



^'7 



-dl(0(a-x)2x) -Y, 5 ^ 



X 



fj{x)fj{a) 



25(x — a) + 2a5'(x — a) — i5(x — a) 



ff(a) 




— d x (6(x — a)2a) 
p > 



5(x — a) 



2(n 



- 1 



S(x — a) 



step 
p = 1 . 



(B9) 



The step leading to the one but last line is a bit dubious for p > 0, since in the second term on 
the preceding line 5(x — a) multiplies a discontinuous /j(x). However, for p > the sum rule 
merely states that J2jfH a ) = 2/p, which can be verified without any distributional 



trickery by contour integration of uG R (a, a; 
converges quadratically, hence also J2j fj( a ) e 



Since fj 



0(j 



Ejfj(a)fj(a + e) 



if p > the sum 
[e i 0) = 2/p. On 
the other hand, if p > and x < a one has J2j fj( x )fj( a ) — 0- 

Returning to g(x), it is straightforward to repeat the calculation for a small positive t, 
leading to 



E^( a+ - x )p( x )fj( x )fj( a ) e " 



,0+ 



' 25{x 
2rT 



rr + 1 
_ 5(x — a) 



a) , |U > 

5 (x — a) , step 

p = 1 



(B10) 



Comparing the outcome with the familiar case y < a, one sees that the presence of the 
^-function has no effect if p > (since the whole weight of the surface term in p falls inside 
[0, a + )). However, if p(x) is continuous at x — a, the factor 9(a — x) halves the weight of 
the QNM sum, as could be expected. The case of a step in p(x) at x = a is intermediate 
between these two (for n~ > 1). 

For reference I also give the "unrestricted" form: 
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Ep(*)mfM=-4 fj{x) { j{a) 
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5(x — a) 



S(x — a) 



fj(x)fj(a) 



Lo'l 



— cL 



fj(x)fj(a) 



10 a 



2-i 



fj{ a )fj( a + e ) 



UJi 



(e|0) 



fl > 



' 2S(x — a) , 

6(x — a) , step 

n~ — 1 

k divergent , 



(Bll) 



p — 1 



However, this last sum rule is only valid up to a finite regular contribution for x > a, so it 
can only be used inside integrals J a dx. 

Finally, I consider an object which is one step more singular, defining r(x,y) = 
J2j^j0(a + — x)9(a + — y)p(x)p(y)fj(x)fj(y), which one for instance encounters upon substi- 
tuting the projection formula for ctj into <j) = djfj. For variation, I present the calculation 
using a test function: 



/ dxdy p{x)p{y)f j (x)f j {y)(j){x,y) 

= ~Y. f dxdyp(y) ^ -f-A-L 



53 / 



+53/ dxdyp(y) 



CO a 



d x (j)(x, y) 



-i£ r + dyp{y)f 3 {a)f 3 {y)<P{a,y) + £ /""dy p (y) lMlM ^ (o>y) 



I] / dxdyp(y) 



fj( x )fj(y) 02 , / \ 

-d x <j>{x,y) 



10 A 



2i(n 



a, a 



(B12) 



r(x,y) 



2ifn~ 



1 



<5(x — a)5(y — a) . 



(B13) 



In the step leading to the third line of (|B12j ) we supposed our test function to satisfy 
0(0, y) = 0; in the one but last line of ( |B12j ), the second and third term vanish on behalf of 
( [BID for < x < y < a. 

Eq. ( |B13| ) leads one to suspect that, for p > 0, J2jUjfj(a) = —2i/p 2 . To verify this, 
integrate [to 2 /(l + u 2 /K 2 )]G R (a, a; u) over the circle \u\ = R. As R — > oo for fixed i^, the 
contour integral vanishes. Now let — > oo; the QNM pole contributions can be evaluated 
for K = oo, and the residues at lo = ±iK can be evaluated in the WKB approximation. 
Combination gives the desired result. 
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With only a trivial variation, one obtains 

2in~ 



■S(x — a)5(y — a) , step 



£ oj j 9(a+ - x)6(a + - y)p(x)p(y)f J (x)f J (y)e-^ 0+ = { n~ + l 

-iS(x — a)S(y — a) , p = 1 . 

(B14) 

However, for p > I find ^2jUjfj(a)e~ lu ' i0 = —2i(n~ + l)/p 2 (the effect of the regulator 
can be found using a WKB calculation of the (^(j^ 1 ) tail of the summand). Apparently, the 
highly singular regulated sum ( |B14|) with its worse convergence no longer allows commuting 
the limit n~ — » oo with the sum over j. 



APPENDIX C: DIAGONAL TWO- VARIABLE QNM EXPANSION 



Let us try to have a more general look at "diagonal" expansions such as (|5.3|) and 
( |5.9| ), which also occur several times in Ref. 0. Thus, suppose one has the one- component 
expansion 



(f>{x,y) = J2 a ifi( x )fj(y) 



(CI) 



for a certain symmetric <fi on (x, y) E [0,af. One sees at once that [p 1 (x)d 2 — 
p~ l (y)dy] <p(x, y) = 0, so that (f> is completely determined by a set of boundary conditions 
which specify a unique solution to a hyperbolic equation. Hence, {fj(x)fj(y)} is grossly 
undercomplete in the space of functions on < x < y < a. 

The existence question for the expansion (|C1[) thus has a negative answer in general. 
However, one has uniqueness under the assumption that J2j a j converges absolutely. For a 
proof, suppose that = 0, i.e., 



(C2) 



Operating on (|C2|) with J a dx p(x) fk(x) yields 



j 

Evaluating ( P3"D at y = a gives 



5 jk - i 



fj(a)fk(a) 



lit 



: fj(o)fk(a) 



while operating on ( |U3|) with —id y \ a - shows that 



.fj{a)fk(a) 

UJj + LO k 



(C3) 



(C4) 



(C5) 
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Furthermore, — p 1 (a )d% 



leads to 



5 jk - i 



fj(a)fk(a) 

UJj + UJ k 



and evaluating ( |C2"D at x = y = a one gets 

2>tf(o) = t>. 



(C6) 



(C7) 



Finally, + @ + z/iQ + z/*(a)(03) reads 

2uJkfk{a)a k = =>- a fc = 



□ 



(C8) 



The technical condition on the summability of aj enables the limit y f a to be taken behind 
J2j to arrive at (|C^)-([U6|). Namely, for jj, > the factor in square brackets is 0(j~ 2 ), and 



this is multiplied at most with ajUi^fjiy) where fj(y) is bounded. Hence, J2j converges 
uniformly with respect to y. If \i — the factor in square brackets is but now the 

prefactor is at most ajUjfj(y) since ( |C6|) is not needed. 

We have not exhaustively examined the intricacies of poorly converging or distributional 
expansions flCT|). However, the above theorem and our experience suggest very strongly that 
the only freedom in the expansion consists in the addition of a term aj = c/uj if /x > 0; 
since in the course of the derivation we also supposed convergence to y) at the point 
x = y = a, by (|B4[ ) even this freedom is absent if /i = (step discontinuity). 
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